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Abstract 

We develop the effective non-Hermitian Hamiltonian approach for open sys¬ 
tems with Neumann boundary conditions. The approach can be used for cal¬ 
culating the scattering matrix and the scattering function in open resonator- 
waveguide systems. In higher than one dimensions the method represents acous¬ 
tic coupled mode theory in which the scattering solution within an open res¬ 
onator is found in the form of expansion over the eigenmodes of the closed 
resonator decoupled from the waveguides. The problem of finding the transmis¬ 
sion spectra is reduced to solving a set of linear equations with a non-Hermitian 
matrix whose anti-Hermitian term accounts for coupling between the resonator 
eigenmodes and the scattering channels of the waveguides. Numerical applica¬ 
tions to acoustic two-, and three-dimensional resonator-waveguide problems are 
considered. 


1. Introduction 

The approach of the effective non-Hermitian Hamiltonian 0,ii have found 
numerous applications in various branches of physics including atomic nuclei 

a s, chaotic billiards SSSSEEl, tight-binding models BIIIH 
Tlijl . potential scattering [^, photonic crystals [l^, etc. The objective 
of the present paper is to revisit the concept of the effective non-Hermitian 
Hamiltonian in application to open resonators with the Neumann boundary 
conditions. The problem of resonant scattering typically involves a resonator 
(which could be an atom, atomic nucleus, quantum dot, microwave or acoustic 
cavity etc) and one, two or more scattering channels that couple the resonator 
to the environment. The mainstream idea is to split the full Hilbert space 
into subspaces: subspace B formed by the eigenfunctions of discrete spectrum 
localized within the scattering center, and subspace C which spans the extended 
eigenfunctions of the scattering channels. Therefore, the exact description of 
open system meets a problem of matching the wave functions of discrete and 
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continuous spectra. In 1957 Livshits [l^ and independently Feshbach in 1958 
[ 2 ^ introduced the idea to project the total Hilbert space onto the discrete 
states of subspace B. Given the Hamilton operator of the whole system as 

+ ^bc + ^cb) (1) 

c 

the projection onto the discrete subsmce leads to the concept of the effective 
non-Hermitian Hamiltonian 

'Heff='HB + 'y^VBC - ^VcB- ( 2 ) 

^ E+-'Hc 


Here 'Hs is the Hamiltonian of the closed system, T-Lc is the Hamiltonian of the 
scattering channel C, Vbc^ Vcb stand for the coupling matrix elements between 
the eigenstates oi'Hs and the eigenstates of the scattering channel C, and E is 
the energy of scattered particle (wave). The term E^ = E+iQ ensures that only 
outgoing waves will be present in the solution after the scattering occurs. As 
a result the effective Hamiltonian Q is a non-Hermitian matrix with complex 
eigenvalues z\ which determine the energies and lifetimes of the resonant states 
as Re{z\), and —2Im{z\) correspondingly [l|, |^. If the propagation band of 
waveguide is infinite then the effective non-Hermitian Hamiltonian takes the 
most simple form widely used in the scattering theory [H, 


111 


He// = HS - fTT ^ VFclT^ 


( 3 ) 


c=i 


where Wc is a column matrix whose elements account for the coupling of each 
individual inner state to the scattering channel C, and the symbol f stands for 
Hermitian transpose. The scattering matrix Sec is then given by the inverse 
oiE-n,ff 1,11 

Scc{E) = Scc-^7riWl, -1 .—Wc (4) 

E-Ueff 

where dec is the Kronecker delta. 

The approach of the effective non-Hermitian Hamiltonian for two dimen¬ 
sional resonator-waveguide systems controlled by the Schroodinger equation 
was previously addressed in Refs. [iii- In particular in Ref. @ the au¬ 
thors derived exact formulas for the coupling matrices Wc for both Dirichlet 
and Neumann Boundary conditions on the boundary of the resonator. What is 
more interesting, however, it was shown Q that only in the case of the Neu¬ 
mann boundary condition the approach is stable with respect to truncation 
of the discrete basis to a finite number of eigenstates thanks to the absolute 
convergence of the spectral sum for the reaction matrix. In the above refer¬ 
ences H, H the resonator-waveguide problem was considered in the context 
of quantum scattering. This imposes restrictions on the applicability of the ef¬ 
fective non-Hermitian Hamiltonian because one normally requires the Dirichlet 
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boundary conditions on the infinitely high ward-wall boundary. Although, there 
are techniques to improve the convergence of the spectral sum [2l|, [22| in the 
Dirichlet case, in general one would resort to the mixed boundary conditions 
applying the Dirichlet boundary conditions on the physical boundaries while the 
Neumann boundary condition is applied on the waveguide-resonator interface 
Q. One the other hand, the Neumann boundary conditions are the boundary 
conditions for the pressure field on a sound hard boundary. That prompts us to 
apply the effective non-Hermitian Hamiltonian formalism to acoustic scattering 
problem. In essence, the proposed approach relying on the spectral properties 
of a closed resonator decoupled from the environment is analogous to the cou¬ 
pled mode theory ^ 241 which is a very popular tool for analyzing resonant 
scattering in optics. Technically, the optical coupled mode theory 2^ 241 rep¬ 
resents a method for finding transmission spectra from a set of linear equations 
with a matrix analogous to Eq.fjH]) in which the diagonal Hermitian term con¬ 
sists of the eigenfrequencies of the resonator, while the second ani-Hermitian 
term accounts for the coupling between the eigenmodes of the resonator with 
the scattering channels. In this paper we will focus on developing acoustic cou¬ 
pled mode theory including applications to discretized systems which render the 
method applicable to open acoustic resonators of arbitrary shape in which the 
eigenfunctions could not be found analytically. 

The article is organized as follows. In Sec. [5] we consider a simple one¬ 
dimensional tight-binding model which is aimed to illustrate our approach to 
derive the effective non-Hermitian Hamiltonian. In Sec. [3]we extend our results 
to 2D case and demonstrate the connection between the discrete model based on 
the finite-difference representation of the Helmholtz equation and the continuous 
model based on the eigenfunctions of the partial differential equation. In Sec. [3] 
we demonstrate an application of the acoustic coupled mode theory for finding 
the transmission spectra and scattering functions in a realistic 3D structure. 
Finally, we conclude in Sec. [5] 


2. The effective non-Hermitian Hamiltonian for one-dimensional sys¬ 
tem with the Neumann boundary conditions 


There are many ^proaches to establish the effective non-Hermitian Hamil¬ 
tonian formalism (USSS 12, 13, 17|. In this paper we adopt a variation of the 
method recently developed for two-particle Bose-Hubbard lattice model [23 . To 
describe the method we begin with the simplest possible one-dimensional tight- 
binding model which consists of one-dimensional resonator coupled to one or two 
half-infinite wires (waveguides). The systems under consideration are sketched 
in Fig. [T]We will see later that this approach can be easily generalized to higher 
dimensions as well as applied to the continuous limit which corresponds to 2D 
and 3D acoustic problems. 
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Figure 1: (Color online) One-dimensional tight-binding chain (resonator) of N sites (red cir¬ 
cles) coupled to one (a) or two (b) semi-infinite waveguides (blue stars). The coupling between 
the resonator and the waveguides (dashed lines) is controlled by hopping matrix element t. 
All other sites are coupled via unit hopping matrix element (solid lines). 


2.1. Tight-binding chain coupled to a single waveguide 

For the model depicted in Fig. [I](a) the ID stationary Schrddinger equation 
takes the following form 


where 


while 



-(2- 

- ^N-l)'^3 = j = ■ • 

.,iV (5) 

f ) = . 

{; 

if j = 0; 
otherwise, 

(6) 

r 

t 

if j = -1; 


= 

0 

if j = V — 1 

(7) 

1 

1 

otherwise. 



One can see that the parameter t in this model controls the coupling between 
the resonator and the waveguide. Meanwhile, because the Neumann BC at site 
N in the discrete case are written as 


ipN — V'Af-i = 0; (8) 

the term is introduced at the end of the tight-binding chain to account 

for the Neumann boundary condition. 

The key idea of our approach is to split the total space into two subspaces. 
The first subspace corresponds to the resonator with j = 0,1,..., iV — 1 while 
the second subspace corresponds to the waveguide with j < 0. Respectively we 
present the solution of the Eq. ([HI) in the following form 


V'(j) 


^bU) 


4 


if j = 0,l,...,iV-l; 

i<o, 


(9) 







where '0 b(j) is the solution within the resonator 

N 

^B{j) = X! XmfpmU), (10) 

m—1 

with tpmU) the eigenfunctions of the resonator. Notice, that here Xm are 
unknown coefficients which are yet to be found for Eq.® to be the solution of 
the discrete Schrodinger equation ([S]). The eigenfunctions are given by 

the following equation. 


^i(j) 




cos(fcm(j+l/2)), 


with the corresponding eigenenergies given by 


7r(m — 1) 
N ^ 


m = 2,3,..., N. 

( 11 ) 


E^n — 2 2cos(fc^). 


( 12 ) 


We note in passing that '(/'m (j) Eq. (fTTll are the eigenvectors of the NxN matrix 
Hamiltonian operator with the Neumann boundary conditions (jS]) at j = 0 and 


j = N-l 
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(13) 


The solution in the waveguide is written as 
outgoing waves 


V'(±)(j) 


1 

y^dTT sin(A:) 


a superposition of ingoing and 

(14) 


with as the incoming and outgoing amplitudes. The waves 'ip±ij) are 
eigenfunctions of the tight-binding infinite wire ~ 1) ~ V'(±)(j + 1) + 

2V'(±)(j) = with the dispersion relation 


E = 2 —2 cos(fc). 


(15) 


Obviously, in this present model the iS-matrix contains only reflection coefficient 
S = r with |r| = 1, i.e., only the phase of the wave changes in the scattering. 

In order to find the equation for coefficients Xrm we evaluate the scalar 
product {ipml'E — E\'ijj)=0 with operator H implicitly defined through Eq. JS]). 
We mention in passing that the bra-ket notation will be adopted thr oug hout the 
paper for vector quantities to comply with the previous works [ij,ll7| . Specif¬ 
ically, I...) stands for a column vector, while (... | = |.. .)^ is the corresponding 
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Hermitian conjugate row vector. 


N-2 


(V’mI'H - E\'ip) = (2 - E)Xm - E V’m(j)'0(j + 1) 

3=0 


N-1 


- E - 1) - IpniiN - l)tlj{N - 1) - tV'm(0)V’(-l) 


= (2 - £;)Xm - tV'm(0)^/’(-l) 


w 


- T, Xr, 


n—1 


N-2 


N-1 


E V'mO')V’n(j + 1) + E “ 1) “ “ 1) 


N ^ N 

= -EXm + E Xn(V’mlHjvlV'n) + E Xn V”™ (0) V'n (0) “ (0)'!/'(-1) 

n=l n—1 


N 

= {Em — E)Xm + E Xn'lprn{0)lpn{0) — ~/^^^'4’ni{0) 

n—1 ^ 


a^+^e-^'^ + a[ 


.(16) 


On the other hand the scalar quantity (^(_)|H — E\'tjj) = 0 could be evaluated 
as 


(V^(_)|7/ - E\^|J) = [^(-2) + - (2 - E)i;{-1)] 


-2 


- E V’(-)(j) [V’(j -1) + +1) - (2 - ^^)V’(j)] 

j^-oo 


477 sin k 


aW+a« 


g-zfc _ 


N 




:e E Xn^/’n(l). 


n=l 


(17) 


The resulting expressions could be presented in a matrix form as 
4-—-e*'=lTL 


Em-E + WLWl 

^ v47rsinfc 

_ t 1 y^ — ik 


\/47r sin k 


lx) 

J-) 




(+) 

L 


— / x/At: sin k 

_, , M,r 1 1 1 -ik (+) 

Air sink ) ^ L J 47r sin fe“L 

( 18 ) 

where |x) is the column vector of coefficients Xm which define the scattering 
function ij^Bij) through Eq. (fTUl) . while Wl is a column vector with its elements 
given by 

= V'm(l)- (19) 

From the second row of Eq. m we have 

+ ty47rsin(fc)w2lx)- 

Substituting this in the first row of Eq. we obtain 


{Ueff - E)\x) = -it\ - sin(fc)VUa 


1 . (+) 




( 20 ) 


where operator Heff could be easily recognized as the effective non-Hermitian 
Hamiltonian 

neff = Em5mn + {l-t^e^'^)WLWl ( 21 ) 
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2.2. Transmission through one dimensional system 

The model considered in the previous subsection was introduced with the 
only goal to illustrate our approach to derive the effective non-Hermitian Hamil¬ 
tonian. Let us now consider a less trivial example, namely, the tight-binding 
model depicted in Fig. [I](b). Repeating the calculations from the previous sub¬ 
section one can easily find that the effective non-Hermitian Hamiltonian now 
takes the following form 

-Heff = EmSmn + {l-t^e^'^ ) ^ WcW^. 

C=L,R 

with 

The equation for the scattering function now reads 


(■^e// - E)\x) = -it)J^ sin(fc) ^ 


( 22 ) 

(23) 

(24) 


C=L,R 


while the reflection amplitudes can be found as 

«c”^ = + \/47r sm{k)Wl;\x) 


(25) 


The iS-matrix is implicitly defined through the following equation connecting 
incoming and outgoing amplitudes 


4 -^ 

d-) 


r t' 
t r' 


4"^ 

d+) 


(26) 


where \ ^ are the right and left outgoing amplitudes, r, r' are the reflection 

coefficients, and t, t' are the transmission coefficients. Combining the above 
equation we find the 5-matrix in the following form is 


Sc,C' = -fc.C' - ityj2s\n{k)Wl; 


Ueff - E 


t\/2sin{k)Wc’ 


(27) 


After some algebra it could be easily found that in the coordinate represen¬ 
tation Hef f takes the following form 

(28) 

where 77 _d is the tight-binding Hamiltonian of the chain of length N with the 
Dirichlet boundary conditions 12, 1^ 
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It is also worth mentioning that the effective non-Hermitian Hamiltonian in the 
Dirichlet representation has the following form [l^ 

=EmSmn-t^e^’^ Y, WcW^c (30) 

C=L,R 


where{IHL}m = V'm(O), {W R}m = ■'^ith V’m(j) as the eigenfunctions 

of the Dirichlet eigenvalue proplem. One can see that Eq. (E51) is fully consistent 
with the earlier findings |12lll3l|. However the present approach is much simpler, 
in particular it is free from evaluation of principal value integrals 


Finally, to illustrate the effective non-Hermitian Hamiltonian approach we 
plot transmission probability T vs. energy E in Fig. [51 To obtain the trans¬ 
mission probability T one can directly apply Eq. (EZl). However, this approach 
would be numerically inefficient because it involves evaluation of the inverse of 
the matrix Ti-eff — E. In practice to speed-up the performance one would first 
use Eq. (l24l) to find the interior wave function jy) and then use Eq. (l25|) for 
finding the reflection amplitudes. The data in Fig. [2] obviously coincide with 
the transmission probability through tight-binding chain with the Dirichlet BC 
according to [l^ . The reason for this is that the resonator-waveguide bound ary 


is artificial and, thus, one has freedom in choosing the boundary conditions 17 [ 
as long as the resulting equations (1551) is consistent with the total Hamiltonian 


of the system ([S]). However, as it will be seen below, the transmission coefficient 
depends on the choice of the boundary conditions at physical boundaries of the 
resonator. 


3. Two-dimensional model 


Before proceeding to the effective non-Hermitian Hamiltonian we would like 
to make some remarks to draw an analogue between the tight-binding model 
and acoustics. The stationary acoustic field satisfies the Helmholtz equation 

V2V’+^V' = 0, (31) 

where Tp is the pressure, c - the speed of sound, and w is the frequency. In what 
follows c is set equal to unity. The condition on the sound hard boundary is 
the Neumann boundary condition which requires the normal derivative of the 
pressure field to be equal to zero on the boundary S 


dn 



(32) 


In this section we consider the resonator shown in Fig. [31 This system can 
be seen as a limiting case of three-dimensional planar duct-cavity-duct acoustic 
system with a rectangular cross-section if the thickness of the resonator along 
the 2 -axis is substantially smaller than the length and the width in the x and 
j/-dimensions. 
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To find the numerical solution the Helmholtz equation (IHTl) could be written 
through the centered second difference scheme as 

1 

— (V'i+ij + i’l-i,] + V'tt+I + V'tt-I - + —V’ti = 0, (33) 

Oq c 

where oq is the step-size in space. The resulting expression (l33)) could be viewed 
as an eigenvalue problem 

= Ell:, 

for the tight-binding ’’Hamiltonian” H with four off-diagonals which carry the 
nearest neighbor hopping matrix element l/og , while the ’’energy” is given by 

E = uj^lc'^. 

Notice that Eq. (IMl) preserves all symmetries of the system in Fig. [3] in the 
discretized coordinates Xj = aoj,yi = a^l. therefore one can easily introduce 
the discretized Neumann boundary conditions (|S]) as 

'4’±Ni,j - = 0, i’i,±Ni - = 0 

on the vertical and horizontal boundaries correspondingly (see Fig. [3]) with Ni 
and Nj being the indexes of the edge sites. At this point we would like to 
notice that our choice of the system and the discretization method are solely 
aimed at simplifying the ensuing algebra. If the configuration of the system 
is more involved then the finite-difference scheme would not be the method of 
choice. We will see, however, that in the continuous limit ag —)> 0 the resulting 
expressions are not dependant on the orientation of the grid. On the other 
hand if a more sophisticated meshing procedure, like hp-adaptive finite-element 
method [i^ , is applied one would use the same technique to derive the effective 
non-Hermitian Hamiltonian under the only conditions that the grid is locally 
regular on each waveguide-resonator interface and preserves the translational 
symmetry in the waveguides. 


3.1. The effective non-Hermitian Hamiltonian for 2D lattice model 

The technique described in the previous section can be easily generalized for 
2D case. Let P be the number of sites across the waveguide and N the number 
of sites along the resonator, then the channel functions are given by 

'^LfpiP 0 = J < 0 (34) 

and 

j> N -l, (35) 

where 4>q(l) is the transversal wavefunction of the waveguide 


</>i(0 = 



0p(O 


[2 (^{p-l)(l + P/2)\ 

- p - j 


p = 2, 3, ... P, 


(36) 
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while the travelling waves and are 


^rHj) = 


ao 


: exp {±ikpaoj ), (j ) = 


Oq 


y^dTT sin(fcpao) y^47r sin(A:pao) 

The corresponding dispersion relationship is given by 

'7r(p- 1) 


" a? 


0 L 


4 — 2 cos{aokp) — 2 cos 


exp (zLikpaoj =F ikpa^N). 

(37) 

(38) 


Using (1551) one can easily check that the waves dSZD obey the normalization 
condition 


Emcux 

01/)= J =-^ J = 6jf. (39) 

Emin 


The scattering domain is a rectangular box with N sites in the direction 
parallel to the transport axis and M sites in the transversal direction. Then the 
eigenenergies of the box with the Neumann boundary conditions are given by 


, „ , 7r(n — 1) \ „ / Trim — 1) 

4 — 2 cos ( - ——- I — 2 cos ' 


N 


M 


(40) 


while the eigen-functions il}n,m{jy) are the products of two factors each given 
by Eq. (fTT|) 

1) = 0«(j)/m(0 (41) 

After evaluating expressions analogous to Eqs. dSl) and dnl) one finds that 
the coupling between m,p mode of the box and q channel in the left waveguide 
is accounted for through the matrix WL^p{n,m) 


WL,p{m,n) = 


(42) 


while for the right waveguide we have 


WR.p{m,n) = ^ V’m(0</'p(0- 

Vao ^ 


The effective non-Hermitian Hamiltonian reads 

^ ^ — pi^pO-a\ 

He// = E„,„<5/C'+E^-^ E Wc,pWl^. 

P=1 “0 C=L,R 

The equation for the interior wave-function \x) is 

p 


i'Heff-E)\x) = -*E 


I sm(kpao) 


E Wc,: 


p=l V C=L,R 


„(+) 

P'^C,p^ 


(43) 


(44) 


(45) 
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and, finally, for reflection amplitudes Qq ^ we have 


(-) J+) , / 47rsin(fcpao) ^^^t 

ao 


“c,p — ®C,p 


(46) 


3.2. Continuous limit 

In the continuous limit uq —>■ 0 the eigenfunctions of the resonator 
i’n,m{x,y) = 'lj^n{x)lpm{y) 

are given by 

' ^ 2 


^i(a^) = \/7Z> tpnix) = ^1 f-cos 


V'i(y) = 


1 


’^miy) = 


cos 


'Km{y-\-Ly{2) 

Ltj 


n > l,x e [0 Lx], 
m>l,y€[-Ly/2 Ly/2^7) 


where we assumed that the a;-axis coincides with the center-line of the waveg¬ 
uides. The corresponding eigenvalues are 


B — 


7rn\ / Trm 

Lx) \ Ly 


The the solutions in the waveguides are 

1 




\/47rfc^ 


e^"'"^^(l>p{y), X <0, 


and 




1 


X > Lx, 


where (jjqiv) is the transversal wavefunction of the waveguide 

'irpiy + d/2) 


^i(y) = \h^ Mv) = 1/3 cos 


The corresponding dispersion relation reads 


, p = 2, 3, ...Q. 


Ep — kp 


/TTp\'^ 

\~d) ■ 


(48) 

(49) 

(50) 

(51) 

(52) 


Taking into account that the step-size oq is found as 


P M iV 


(53) 


the continuous limit from Eqs. (El, dSl, and dn to Eqs. dUD, dH, and dSOl) 
is reached as 


i’alix,y)= lim 

ao->0 


V'c|'p(j,0 ] , i’n,m{x,y) = lim 


1 


ao—>-0 y ^yCLo 


'0n,m (ji 0 ) ; 


(54) 
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Thus, the continuous limit of coupling matrices Wc,q{m,p) could be found as 


d/2 

WL,pin,m) = / dy^prn{y)(t>piy) 

-d/2 


(55) 


while for the right waveguide we have 


d/2 

WR^p{n,m) =tl:n{Lx) / dyiprn{y)(l)p{y), (56) 

-d/2 

then the effective non-Hermitian Hamiltonian reads 

CXD 

^e// = ~ ^C,pWQ p. (57) 

p=l C=L,R 

The equation for the scattering function takes the following form 

oo r, 

(He//-i?)|x) = -*EVV ^ (58) 

p=l ’ C=L,R 

while the reflection amplitudes are given by 

W^Jx) (59) 

The interior wave function 'ipi,(x,y) could be found as 

1pb{x,y) = ^ Xn,m'tpm,ni^^y) (60) 

n,m—l 

The above presented approach is valid not only for the system shown in Fig. 
[3l In general the resonator eigenfunctions could not be found through separation 
of variables as in Eq. (j35p . However, on any zero curvature waveguide-resonator 
interface of length d one can introduce a local set of coordinates with the y-axis 
directed along the interface while the x-axis is directed along the waveguide. 
Then the same calculations give us the the following expression for the coupling 
matrix 

d/2 

Wp{m)= / dy'ilJm{y,x)(j}p{y), (61) 

-d/2 

where we have a single index m to enumerate the eigenfunctions of the resonator. 
The rest of the formulas (1551 I59I60I) remain correct with the only difference 
that index C could run over arbitrary number of waveguides attached to the 
resonator. This resulting expressions are consistent with the earlier findings in 

ref. 0. 
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3.3. Numerical results 

To perform a numerical test we calculated the transmission coefficient T 
using both discrete (Hil) and continuous (1551) approaches. The geometric pa¬ 
rameters of the system were chosen as d = 1, Ly = 1.8, = 1. In the discrete 

model we took N = 20, Q = 20, P = 36 which corresponds to Oq = 0.05. In the 
continuous model where we have an infinite number of modes in the resonator 
as well as in the waveguide Eq. (I58p was truncated to a finite number of un¬ 
knowns. To be consistent with our discrete model we took 20 modes in each 
waveguide while the indices m and p in Eq. (15811 run as m = 1,2,... ,M and 
n = 1 , 2 ,... , N. Thus, we guaranteed that in both approaches we have the same 
number of degrees of freedom. For our system the coupling matrices WL,p{n, in) 
and Wft,p{n,m) could be found analytically through Eqs. dSS]) and (1551) . The 
resulting expressions are 


WL,pin,m) = 


12-60 


^VL,p{n,m) — 


12-60 


(-irv,, 


(62) 


where 



Vp,m — ^ 


(63) 

The numerical data are shown in Fig. |3]where one can see a very good agreement 
between two methods. This, in accordance with ref. Q, proves that the effective 
non-Hermitian Hamiltonian is stable with respect to truncation of the number 
of modes. Finally, in Fig. [5]we present two typical configurations of the real part 
of the pressure field with the resonator excited through both first and second 
scattering channels. 


4. Transmission through cylindrical acoustic resonator 

In this section we consider sound transmission in three dimensional system, 
namely a cylindrical resonator coupled with two cylindrical waveguides as shown 
in Fig. [HI We restrict ourself to the coaxial case where the angular momentum 
-i-^ is a constant of motion cj) being the azimuthal variable. Therefore the 
Hilbert space of the total system is a direct sum of the subspaces specified by the 
angular momentum. Below we present the effective non-Hermitian Hamiltonian 
approach in case of zero angular momentum. The generalization to non-zero 
angular momentum is straightforward. 


, m = I 
, m = I 
, m > I 

, m > 1 
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4-1. The effective non-Hermitian Hamiltonian 

The eigenfrequencies of the closed cylindrical resonator with the Neuman 
boundary conditions are given by 


bJ 


2 

mn 


= C 




■nvf 

IX 


(64) 


where R and are the radius and length of the cylindrical resonator, and 
is the m-th root of the equation = di(/im) = 0,m = 1 , 2 ,... for the 

derivative of the first order Bessel function. The corresponding eigenfunctions 
are 



(65) 


Similar to the 2D acoustic case considered in the previous section in the con¬ 
tinuous limit we obtain the effective Hamiltonian in the same form as Eq. dSll) 
with the coupling matrices given by 


l)z 


WL,p{'^,m) = 


2 - 6 ° 




^-T^(-l)”-V,,^, (66) 


with 


a 27r 


Vp^m = —t / / drd(j)rJo{^ijnr/R)Joinpr/a), (67) 
T:aRJo{tim)Jo{t^p) J J 
0 0 

where a is the radius of the cylindrical waveguides. Integration in Eq. ([57|l can 
be done analitically [27 1 


^ 2a fjjjYi R) 

inmaf - inpR)^ Joinm.) 

Finally, the dispersion relation for the scattering channels reads 

X = cHXp/a^ + kl),p = 1,2,.... 


( 68 ) 


(69) 


The numerical procedure for finding the reflection amplitudes OpJ and the in¬ 
terior scattering function is identical to the 2D case considered in the previous 
section Eqs. dSB, dsn (IMI), and (EOl) 


4 . 2 . Mixed representation 

Another interesting opportunity for separable systems is to use the dis¬ 
cretized representation in one coordinate while the continuous representation 
is used in the other coordinates. Let us write the Helmholtz equation (1511) in 
the following form 

ipj - 2ijjj + + i/jj+i ] -I- Etjjj = 0, (70) 


1 (\d^ Id' 

Oq \ dr"^ r dr 
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where we used the continuous representation for the radial and azimuthal co¬ 
ordinates and the discrete representation for the axial coordinate z. The idea 
of the approach is sketched in Fig. [71 The resulting expressions are the same 
as those found in lO for both scattering function Eq. (1451) and reflection am¬ 
plitudes Eq. (HSI) with the only difference that the coupling matrices are now 
given by 


WL,p{m,n) = WR^p{m,n) = 

-v/oo Jao 


p,m 5 


(71) 


where Vp^m are defined through Eq. (1551) . The eigenfrequincies of the resonator 
are now found as 


Wmn = 


2-2cos(fc„ao) 

' o 


a, 

while the dispersion relation for the channel function is 

4 

a? a. 


,fc„ = = (72) 


2 2 
(jJ = C 


Mp 2 — 2 cos(fcpao) 


0 


,P= 1,2,.... 


(73) 


4 . 3 . Numerical results 

Numerically computed scattering function are shown in Fig. [51 Here, to 
perform the numerical test we took 10 modes in each waveguide while in the 
resonator 20 modes were taken in each radial and axial variables. The same 
numbers were chosen in the mixed representation. In Fig. j^lwe demonstrate the 
transmission spectra of the cylindrical resonators with two differing aspect ratios 
a/R. One can see a good agreement between continual and mixed approaches 
even for relatively small number of slices TV = 20 (see Fig. O. 


5. Summary 


In this paper we developed the approach of the effective non-Hermitian 
Hamiltonian to resonator-waveguide systems with the Neumann boundary con¬ 
ditions. To derive the effective non-Hermitian Hamiltonian we used an alter¬ 


native approach which dates back to the seminal paper by Fano 28| and was 
recently adapted for the discrete Schrodinger equation [ 2 ^ . In contrast with the 
standart Feshbash projection technique [20| our approach does not face a prob¬ 
lem of principal value integration. It was shown that the method of the effective 
non-Hermitian Hamiltonian can be used to calculate the ^-matrix and the scat¬ 
tering function of sound-hard boundary acoustic scattering problem. Thanks to 
the projection of the total Hilbert space onto the inner states of the resonator 
the method is free of the discontinuity problem typical for the mode-matching 
techniques [ii, [IS im, m 113, in m, 


36l |. It was shown that in the continuous 
case the effective Hamiltonian can be truncated to a reasonable small num¬ 
ber of eigenmodes of the closed system whose eigenfrequencies are concentrated 
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around the frequency of the incident wave. In that sense the method of effective 
non-Hermitian Hamiltonian is analogous to the coupled-mode theory [^, [ 2 ^ 
which however neglects the radiation shifts caused by a finite lower band edge 
and dispersive properties of the waveguides. Therefore we can refer to the ap¬ 
proach of the effective non-Hermitian Hamiltonian as an advanced form of the 
coupled mode theory. In this paper the effective non-Hermitian Hamiltonian 
was formulated in both continuous and discrete forms. The former is shown 
to be consistent with the earlier findings in ref. @ while the latter seems to 
be more feasible for numerical implementation along with finite-difference and 
finite-element discretization schemes. In fact, the developed approach represents 
a tool for modelling systems with o pen boundaries. Alike to the recently devel¬ 


oped wave finite-element methods |37l . l38j it explicitly utilizes the translation 


invariance in the waveguide to formulate numerically exact transparent bound¬ 
ary conditions [s^. It should be noted, however, that the resulting equations 
(m, (03), and (1581) contain dense rather then sparse matrices typical for the 
lattice methods. This to a certain degree degrades the numerical performance. 
Nevertheless, thanks to a clear physical picture of the wave transmission one 
may think of developing techniques based on truncation of the inner space to a 
small number of modes which are relevant to a particular case of the scattering 
problem. This represents a goal for the future studies. Finally, the effective 
Hamiltonian has an advantage in physical interpretation of it complex eigenval¬ 
ues z. The real parts Er = Re{z) define the positions of resonances with the 
resonance widths defined by —2Im(z) [^. That advantage becomes important 
in study of wave trapping [i^, EH, EH] by tracing the imaginary parts of complex 
eigenvalues of the effective non-Hermitian Hamiltonian. 
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Figure 2: (Color online)Transmission spectra of tight-binding chains of 5 (a) and 7 (b) sites, 
t = 0.4 
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Figure 3: (Color online) 2D acoustic resonator coupled two waveguides. 
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Figure 4: (Color online) Transmission spectra of 2D acoustic resonator with the incoming 
wave in the first (a) and second (b) channels; discrete model - solid blue line, continuous 
model - dashed red line. 
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Figure 5: Scattering function at lj = 37r/2 with the incoming wave in the first (a) and second 
(b) channels. 



Figure 6: Cylindric resonator of radius R and length L connected to coaxial waveguides with 
radius a. 
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Figure 7: Cylindric resonator discretized in the axial coordinate. 



Figure 8: Scattering function (real part) in a cylindrical resonator coupled with two coaxial 
waveguides. The wave is incident through the fist channel. L = 3.5, a; = 6.6845. (a) R = 1.5a 
and (b) for R = 2a. 
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Figure 9: (Color online) Transmission spectra of cylindric acoustic resonator with two differing 
aspect ratios: (a) Kja = 1.5 and (b) Kja = 2. Continuous representation - solid blue line, 
mixed representation - dashed red line. 
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